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We report quantum oscillations in the underdoped high temperature superconductor 
YBa2Cu306+a; over a wide range in magnetic field 28 < fioH < 85 T corresponding to ~ 12 oscilla- 
tions, enabling the Fermi surface topology to be mapped to high resolution. As earlier reported by 
Sebastian et al. [Nature 454, 200 (2008)], we find a Fermi surface comprising multiple pockets, as 
revealed by the additional distinct quantum oscillation frequencies and harmonics reported in this 
work. We find the originally reported broad low frequency Fourier peak at « 535 T to be clearly 
resolved into three separate peaks at ~ 460 T, ~ 532 T and w 602 T, in reasonable agreement with 
the reported frequencies of Audouard et al. [Phys. Rev. Lett 103, 157003 (2009)]. However, our 
increased resolution and angle-resolved measurements identify these frequencies to originate from 
two similarly-sized pockets with greatly contrasting degrees of interlayer corrugation. The spectrally 
dominant frequency originates from a pocket (denoted a) that is almost ideally two-dimensional in 
form (exhibiting negligible interlayer corrugation). In contrast, the newly resolved weaker adjacent 
spectral features originate from a deeply corrugated pocket (denoted 7). On comparison with band 
structure, the d-wave symmetry of the interlayer dispersion locates the minimally corrugated a 
pocket at the 'nodal' point k no dai = (tt/2,tt/2), and the significantly corrugated 7 pocket at the 
'antinodal' point k ant i n odai = (tt, 0) within the Brillouin zone. The differently corrugated pockets at 
different locations indicate creation by translational symmetry breaking — a spin density wave is 
suggested from the strong suppression of Zeeman splitting for the spectrally dominant pocket, addi- 
tional evidence for which is provided from the harmonics we resolve in the present experiments. In 
a broken-translational symmetry scenario, symmetry points to the nodal (a) pocket corresponding 
to holes, with the weaker antinodal (7) pocket corresponding to electrons — likely responsible for 
the negative Hall coefficient reported by LeBoeuf et al. [Nature 450, 533 (2007)]. Given the simi- 
larity in a and 7 pocket volumes, their opposite carrier type and the previous report of a diverging 
effective mass in Sebastian et al. [Proc. Nat. Am. Soc. 107, 6175 (2010)], we discuss the possibility 
of a secondary Fermi surface instability at low dopings of the excitonic insulator type, associated 
with the metal-insulator quantum critical point. Its potential involvement in the enhancement of 
superconducting transition temperatures is also discussed. 

PACS numbers: 74.25. Jb, 74.72.-h, 74.72. Gh, 71.30.+h, 71.35.-y, 71.18.+y 



I. INTRODUCTION 



Different bosonic physics from tightly bound pairs of fermions may arise depending on whether the binding takes 
place between like particles or between particles and holes [1]. Such bosonic physics can potentially be involved in 
the physics of high temperature superconductivity. On one hand, strongly interacting pairs of like particles with zero 
momentum have been proposed to constitute the fabric of unconventional superconductors 0, H[, while on the other, 
electron-hole pairs with finite momentum could condense into a competing state with a superlattice, such as a spin 
density wave. Another realisation of such particle-hole pairing could occur in the limit of strong Coulomb coupling, 
where electron and hole pockets of identical size become susceptible to an excitonic insulator instability [4j-|6|. In 
this study, we use high resolution quantum oscillation measurements to show that the underdoped high temperature 
superconductor YBa2Cu306+2; Q contains compensated electron and hole pockets, potentially predisposing it to 
excitonic electron-hole pair condensation. 

Quantum oscillations have been measured on underdoped YBa2Cu30_6+x using in-plane transport [H out-of- 
plane transport [101 ]. torque llMl4|. and contactless conductivity 15Ml7t| experiments. In this work, the power of 
quantum oscillations is enhanced by the technique we employ that utilises interlayer dispersion measurements to 
enable a momentum-space identification of individual Fermi surface sections— details of which are provided in Section 
II. Quantum oscillations are measured on high quality detwinned single crystals of YBa2Cu306+ x (x=0.54, 0.56) Q 
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as a function of magnetic field and of angle. Measurements are made down towlK with the sample immersed in 4 He 
medium, and the contactless conductivity technique used to obtain a high value of signal-to-noise ratio. Out-of-plane 
rotation measurements are performed in pulsed magnetic fields up to 65 T, and two-axis (in-plane and out-of-plane) 
rotation measurements are performed in DC fields up to 45 T. Additional experimental and analysis details are 
provided in Appendix A. 

Our angle-resolved measurements show that the dominant series of low frequency oscillations 8] arises not from one, 
but from two pockets of similar size, the distinct topologies of which point to different positions in the Brillouin zone 
on comparison with the interlayer dispersion, which has d-wave symmetry in underdoped Yl^CuaOe+x 18, ll| (see 
Fig. [Tk). On comparison with the interlayer hopping integrals in YBa2Cu306+a; [HI, EH, we show that the prominent 
frequency corresponding to the a pocket with minimal interlayer corrugation is likely to represent a section of Fermi 
surface at the (±7r/2, ±7r/2) nodal location in the Brillouin zone where holes are expected to nucleate on doping. In 
contrast, the more recently measured satellite frequencies of smaller amplitude correspond to a strongly corrugated 
pocket (7) close to the antinodal location in the Brillouin zone, suggestive of electron carriers. 
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FIG. 1: a, A schematic of the small hole (blue) pocket a and small electron (red) pocket 7, illustrating greatly contrasting 
degrees of corrugation and approximate locations within the sector of the Brillouin zone bounded by (0,0), (7r/2,7r/2) and 
(•7r,0), shown for — 2tt < k z < 2n. Also shown in semi-transparent fashion is the larger hole orbit /3 (corresponding to a 
separate Fermi surface section or potentially a magnetic breakdown orbit [I3|)- The degree of corrugation is determined from 
angle-dependent Fourier transforms shown in Figs.0and frequencies shown in Fig. 3^,. On using the Onsager relation, the three 
pockets collectively yield an effective hole doping of p = If .7 ± 0.2 %. b, Estimates of the ratio of the exciton binding energy 
to kinetic energy Ryd* / E\^ c and effective exciton separation r s as a function of oxygen concentration x, using the expressions 
in Appendix B and effective masses from Refs. [la. JTrj]. with arrows indicating the appropriate axes. Also shown are the 
corresponding values of T c (renormalized) from Ref. [7J and different shadings representing the insulating and metallic regimes 
at low T. 



II. RESULTS 



Interlayer corrugation 

High experimental resolution of the quantum oscillation frequency is key to resolving and characterising the two 
adjacent peaks flanking the larger central peak (shown in Fig. The weaker spectral features at 460 ± 2 T and 
603 ± 2 T on either side of the central frequency at 533 ± 2 T are distinguished on account of their separation 
being well above the frequency resolution limit of our current experiment (AFii m w [I/A(l/^o-£0] ~ 42 T). The three 
distinct peaks in the Fourier transform (shown in Fig. [2]) are observed both without (black curve) and with (red curve) 
a Hann window to reduce diffraction artefacts. 

Evidence for the three resolved low frequency oscillations originating from similar volume electron and hole pockets 
with different degrees of corrugation (shown schematically in Fig. [1^) in addition to the previously resolved high 



frequency oscillation 12j is presented in Fig. [31 which shows quantum oscillations and Fourier transforms measured at 



different angles 9 between the crystalline c axis and H, and in Fig.[2]where the separately resolved quantum oscillation 
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FIG. 2: a Quantum oscillations made using the contactless conductivity method on a sample of YBa2Cu3C>6.54 over a broad 
range of magnetic fields 28 < ^oH < 85 T at « f.5 K, corresponding to « 12 oscillations of the prominent a frequency. 
The experimental data below and above 65 T are from measurements made using the same sample on the same probe but in 
two different pulsed magnets, b A Fourier transform of the quantum oscillations both without (black curve) and with (red 
curve) a Hann window applied. The inset shows a comparison of the separation in frequency F 7 . nec k — F a or F a — F 7i b c ii y 
with the frequency resolution limit AFn m « [l/A(l//io-ff)] (left axis) and with the equivalent number of a oscillation periods 
n a (bottom axis)— the frequency separation exceeds the frequency resolution for n a > 8, enabling all three frequencies to be 
clearly resolved in Fourier transforms and fits of the quantum oscillations. 



frequencies (F a , F@, i^neck and -F 7 ,bciiy) are plotted versus 9. The subtle beat pattern modulating the amplitude of 
the measured oscillations (Fig. [3^, upper panel of upper inset) and the multiple low frequency peaks in the Fourier 
transform suggest two different Fermi surface sections yielding closely spaced low frequencies. To separate a and 7 
pocket contributions, we use three independent analysis methods and find them to yield mutually consistent results 
for a significant difference in corrugation between the a and 7 pockets. 

Given the layered character of the cuprate family of high temperature superconductors, the quasi-two dimensional 
nature of the Fermi surface needs to be factored into the quantum oscillation analysis. Specifically, for quasi-two 
dimensional materials, interlayer corrugation introduces phase smearing that leads to a distinctive beat pattern which 
has been well characterised [20[. The location and separation in magnetic field of the amplitude zeroes (nodes) 
constituting the beat pattern are a topologically constrained function of the depth of corrugation, and may be used 
to extract the size of the corrugation. For this purpose, we use an expression that explicitly treats the form of phase 
smearing in layered systems f2lj ] (Eqn. [1]) rather than the standard Lifshitz-Kosevich expression [22[ which does not 
capture this form of interference (see Appendix G). 

We perform an unconstrained fit (red line) of two corrugated Fermi surface sections to the measured oscillations 
in the 65 T magnet (black line) for a sample of oxygen concentration x — 0.56 at 9 = (Fig. [3^, upper inset). 
Oscillations are fit to the expression 

A/ = £ A/ 4 , i? T i? D J (^^-) cos (1) 
' Mo-" COS 8' Mo-n 

i— a, 7 



incorporating a summation over independent a and 7 corrugated cylinders 
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23| . The spin splitting factor R s 
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FIG. 3: a, Fourier transform before (purple) and after (magenta) subtracting the a component of a fit to Eqn. \T\ (green 
line in the upper inset) to the quantum oscillations (black line in the upper inset) measured in a sample of Yl^CuaOe+z 
with x = 0.56 at T « 1.5 K using the contactless conductivity technique (see Appendix A). A polynomial is used to fit the 
background and the waveform is modulated by a Hann window prior to Fourier transformation. The prominent oscillations 
and harmonics are indicated. Extremal Fermi surface cross-sections given by Au = (2ne/h)F are obtained from the observed 
frequencies F a , F 7 , n eck and F 7l belly and harmonics 2F a , 2i ? 7incc k and 2f 7i bcii y as indicated; here H — |H|. The red line in the 
upper inset is an unconstrained fit to Eqn. [T] where the green and blue lines are the a and 7 components respectively. The 
purple curve (upper inset, lower panel) shows the result of subtracting the a component (green line) of the unconstrained fit 
to Eqn. \T\ from the raw data, revealing a distinctive beat pattern in the lower panel of the inset. The vertical dotted lines 
represent the fields at which nodes occur for n = 2, 3, 4 from Eqn. [3] using the fit value of A_F 7i o. The lower inset shows the 
temperature-dependences of the F a , F 7 , llcc k and F 7 ,bciiy amplitudes (the latter two renormalized to equalize the amplitudes 
at the highest temperature), indicating similar effective masses m* for the a and 7 pockets. A fit to A/ = aoRr [22| yields 
771^=1.56 ± 0.05 m e , m* neck = m- 7l boiiy — 1-6 ± 0.2 m e (where m c is the free electron mass), b & c, Examples of Fourier 
transforms of oscillations measured in a sample with x — 0.54 at T ~ 1.5 K and different angles 9 (data in Figs. 2p and c) 
before and after subtracting the fitted a component (see Fig. [4b ) as described in the text (expanded to show the finer high 
frequency structure)— further Fourier transforms before and after subtraction are shown in the Appendix D. Note that the 
higher frequency Fp (when observed) is 3.14 ± 0.05 times larger than F a , inconsistent with a harmonic of F a [T2I. [l5|. 



does not appear in this expression since it has been shown to be close to unity in previous measurements (ref. 14| and 
Appendix C). Here, A/^o is the amplitude, Rt = (2tt 2 k-Qm*T / efi^H) / smh(2ir 2 kRm*T / eh^!!) is a temperature 
damping factor where m* is the quasiparticle effective mass obtained from the fit shown in Fig. [31 Ru = exp (^7^-) is 
a Dingle damping factor where Tj is a damping constant [22j . Ai 7 ^ is the depth of corrugation of sheet 'z' at angle 
9, while the Bessel function J captures the interference due to phase smearing, with characteristic 'neck' and 'belly' 
extremal frequencies 

Fi, nock, belly COS 8 = F. t ± AF^g (2) 

for each section n performing an unconstrained fit to Eqn. [T] we obtain corrugations AF a .o 11 T and 

A-F 7j o = 72 ±4 T, while F a .o = 531 ±2 T and F 7i q — 532 ±2 T are the same within experimental uncertainty. These 
parameters correspond to three frequencies F 7incc k = 460 ± 4 T, F a fl — 531 ± 2 T (the splitting being unresolvably 
small) and -F 7 ,beiiy = 604 ± 4 T in close agreement with those resolved in the Fourier transform in Fig. [5J Similar 
frequencies were recently reported in Ref. [13j . in which case the experimental uncertainty was larger due to the 
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reduced range in 1/H (see Fig. [2}d inset). The analysis presented here is performed on data spanning between w 9 
and 12 measured oscillations. 

The monotonic field dependence of the a component (green line) identified in the present experiment implies that 
it can be reliably subtracted from the raw data (on performing a fit to Eqn. [I} , enabling the residual beat pattern of 
the 7 oscillations (purple line) to be extracted (lower panel of the upper inset to Fig. [3^). Distinct nodes are observed 
in the residual 7 oscillations. We turn to the expected locations of zeros in the amplitude due to the interference 
pattern to confirm that they coincide with the observed nodes. The values of magnetic field at which amplitude zeros 
are expected due to interference are given by 



HoH n ,i cosfl = 



An - 



(3) 



Eqn. [T] [21| can be expanded as the interference between neck and belly frequencies Fj — AFj $ and Fj + AFi g 
respectively. The neck and belly frequencies interfere destructively whenever they have a relative phase difference 
equal to an odd multiple of 7r, which occurs at values of the field H n given by Equation [3] 

Using the fit value of AF 7j o, the expected nodes for n = 2, 3, 4 are indicated by the vertical red lines in the lower 
panel of the upper inset to Fig. [3k,. We find these to coincide very well with the observed nodes in the residual 
7 oscillations, establishing that the frequencies F 7jncc k « 460 T and F 7 ,bciiy ~ 604 T correspond to the neck and 
belly extremal frequencies of a single corrugated 7 cylinder, distinct from the a section. Indeed, these frequencies 
correspond to those that appear in the extracted fourier transform in Fig. [3] (depicted before and after subtraction of 
the fitted a oscillation component) and in Fig. [3J 
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FIG. 4: a, The product Fcos9 (for both fundamental and harmonic - solid symbols for unsubtracted data and hollow symbols 
for subtracted data) obtained from the Fourier analysis shown in Fig. [3]d and c and Appendix D plotted versus 9. F 7inG ck and 
F 7i beiiy are obtained by Fourier analysis after subtracting the fitted a component (using Eqn. [TJ. For clarity, the data have 
been symmetrized with respect to ±6: both the swept 9 data (shown in Fig. |6k) and measurements made over extensive range 
of 9 in Ref. [3] indicate the oscillations to be symmetric in ±9, as expected for a detwinned single-phase orthorhombic crystal. 
Lines indicate the expected 9 dependences according to Eqns. [2] and HI using the fitted parameters in Fig. [3Jl Consistency with 
Eqn. [4] indicates a deeply corrugated 7 section of Fermi surface [2CJ, |24| . The harmonics are obtained directly from peaks in 
Fourier transforms of the raw data— consistency between fundamentals and harmonics confirms our finding of similarly-sized 
but distinctly different a and 7 Fermi surface sections, b, Examples of quantum oscillations measured in YBa2Cu306+a; with 
x — 0.54 at different angles 9, as indicated (where <f> w 0) together with the fitted a component (green lines, using Eqn [TJ 
as described in the text, c, The quantum oscillations measured in YBa2Cu306+a: with x = 0.54 at different angles 9 after 
subtracting the green fits in (b), yielding residuals dominated by the 7 pocket oscillations. Curves are offset according to 9 
(right axis). Also plotted (red lines) are the predicted values of H n from Eqns. [3] and [4] using the fit value of AF 7 ,o- 



In Figs. [5}d,c and [H we extend this analysis to a second sample of oxygen concentration x = 0.54, for which 
oscillations are measured at several different angles 6. Fits to Eqn. Q] similar to those in Fig. [3^ are made to the data 
in Fig. gb as a function of angle (here AF Q , ^ 10 T, AF 7 , = 74 ± 4T, F a ,o ~ 535 ± 2 T and F 7 , w 533 ± 2 T). 
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Fig. HJi shows the residual 7 oscillations on subtracting the fitted a component (green lines in Fig. [4]o) , and Fig. [4^, 
shows the extracted peaks in the Fourier transform (shown in Fig. [3^, and b and Appendix D) of the subtracted data 
(hollow symbols) and unsubtracted data (solid symbols). The angular dependent amplitude nodes observed in the 
residual 7 oscillations in Fig. 0t and frequencies in Fig. 0^ are compared with those expected for a single corrugated 
cylinder, by using Yamaji's expression for an angular dependent frequency difference between the 'neck' and 'belly' 



of a single corrugated cylinder [24j: 



AFi, e = AFi, o Jo(ft||Ctan0). (4) 

On substituting the fit value of AF 7i o into Eqn. @J the predicted nodes for a single corrugated cylinder in Eqn. [3] are 
shown by lines in Fig. Hfc, and predicted neck and belly frequencies are shown by red lines in Fig. [4^,. We see that 
excellent agreement with the data is obtained in both Figs.|4ji and c independent of any fitting parameters, confirming 
a single corrugated 7 Fermi surface section (distinct from a). 

The small value of AF a implies that the putative nodes for the a Fermi surface section occur at fields too low for 
interference between neck and belly frequencies to be observed in the form of nodes. We can now understand why 
the a cylinder yields a spectrally dominant peak at F a fl ~ 532 T in Fig. [3^— no significant smearing of the quantum 
oscillation phase occurs on averaging over k z for this section (see diagram in Fig. QJi). The larger corrugation of 
the 7 section, by contrast, implies that significant phase smearing occurs on averaging over k z for this section, with 
phase coherence being achieved only at the extrema -F 7j neck and i^^eiiy - accounting for the spectrally weak satellite 
frequencies flanking the central peak in Fig. [3^. 

For a final consistency check, we turn to the harmonics 2F ai 2i 7 ! 7jncc i t and 2i 7 7j boiiy, which are observed for both 
samples (x = 0.54 and x = 0.56) in Figs|3k, b and c (plotted versus 9 in Fig.[4jt together with the Yamaji prediction). 
Because the harmonics appear in the Fourier transforms of the raw data (Figs. [3^, b and c and Appendix D), they 
provide a second estimate of each frequency (Fig. [3^,) that is independent of any subtraction of the predominant 
a frequency and therefore independent of the previous fits. A further advantage of harmonic detection is that the 
doubling of the frequency yields a twofold increase in the confidence to which the separate Fermi surface sections 
can be resolved. The absolute frequency resolution obtained from the interval in magnetic field is therefore A.Fn m rs 
[l/A(l//io-ff)]/2 ~ 26 T. The lack of a peak splitting or deviation of F a and its harmonic from F a = F a j cos9 over 
the entire angular range in Fig. implies an irresolvably small depth of corrugation AF a .o < 13 T, while the values 
of AF 1: g — (2F 7i bciiy — 2i 7 ' 7i ncck)/4 coincide well with the expected values from the fundamental analysis (shown 
in Fig. @Ji)- The observation of prominent harmonics for both the a and 7 pockets also enables us to estimate the 
relative size of the damping term for each of the electron and hole pockets, which we find to be similar in magnitude: 
r Q ~ T 7 ~ 10 2 T, corresponding to a mean free path of ~ 200 A. 

We therefore find the nodal analysis of the oscillations, Fourier transform analysis, Yamaji fits and second har- 
monic detection all to yield internally consistent results corresponding to two similarly-sized a and 7 sections with 
significantly different degrees of corrugation. Oscillations from two pockets with the same degree of corrugation, as 
expected for bilayer split pocket areas in the underdoped cuprates [13j , on the other hand, would not be consistent 
with our data (see Appendix G). 



Correspondence with equal and opposite electron and hole pockets 



In the simplest model, our finding of Fermi surface sections with significantly different degrees of corrugation 
locates them at different locations in the Brillouin zone. We compare the different measured interlayer hoppings with 
calculated values in a tight binding model of underdoped YBa2Cu306+iE [25| to locate the different pockets. As shown 
in Refs. 13, III the low-energy in-layer dispersion [k = (k x , k y )] of the Cu d x 2_ y 2 — like band is: 



e (k) = -At 



u (k) + 2r 



v(kY 



1 - 2ru (k) 



(5) 



where u(k) = ^(cosk x + cosk y ), v(k) = ^(cosk x — cosk y ), t is the nearest neighbour hopping integral, and r is a 



'range' parameter due to in-layer hopping via the so-called axial orbital 
2p z . For small r, Eqn. [5] reduces to the familiar tight-binding form: 



a hybrid of mostly Cu 4s and apical oxygen 



e (k) = const. — 2t (cos k x + cos k y ) + At' cos k x cos k y 
— 2t" (cos 2k x + cos 2k y ) 



(6) 
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with t' = rt, and t" = hrt. Hopping in the z-direction occurs via the axial orbital and therefore depends on (k x ,k y ) 
as i^2ru(k) ' (^ e- ) ^ nas ^ ne same symmetry as the superconducting gap, vanishing along the nodal lines and reaching 
maxima at the (±7r, 0) and (0, ±7r) points. If the layers are stacked on top of each other, the fc z -dispersion is simply 
included in Eqn.[5]by the substitution: r — > r + (t±/t) cos ck z . For YBa2Cu306+z (x ~ 0.5), t « 400 mcV, tj_ causes 
r to vary between w 0.28 and 0.32 and the corresponding Fermi surface is shown in Fig. [5^. Here the line thickness is 
the size of the k z dispersion yielded by the perturbation in r. The different interlayer hoppings at different locations 
in the Brillouin zone reveals that small pockets at each of these locations created by Fermi surface reconstruction will 
have different corrugations. Accordingly, we proceed to examine the difference in corrugation expected for each of the 
multiple pockets throughout the Brillouin zone yielded for a translational symmetry breaking order parameter. 
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FIG. 5: a, Schematic unreconstructed Fermi surface depicting the translational vectors Q s ,±, the line thickness (obtained by 
a perturbation in r as described in the text) represents the inter-layer dispersion k z responsible for corrugation and is seen 
to be greater in the antinodal (±7r,0), (0,±7r) compared to the nodal (±7r/2,±7r/2) region, b, Representative reconstructed 
Fermi surface according to models used in Refs. fl2j, [26[-[28( (described in the text and Appendices F, I). Blue and red lines 
correspond to the nodal hole pockets and antinodal electron pocket respectively. The line thickness (obtained by a perturbation 
in r) represents the inter-layer dispersion k z and is significantly bigger for the antinodal compared to the nodal pocket, c, 
Reconstructed Fermi surface consisting of only a single type of pocket (/3) after increasing the coupling A s (neglecting A c ). 
Antiferromagnetic bilayer coupling 18] (shown in Fig.QIJ]of Appendix F) would yield two slightly different variants of the Fermi 
surfaces in b and c, but with very similar pocket areas, d, Schematic dispersion prior to Fermi surface reconstruction depicted 
in one-dimension, e, Schematic dispersion of the three pockets, f, Schematic showing the opening of a gap between the a and 
7 pockets, as expected to result from an excitonic insulator instability. Here, we neglect the charge potential A c , the effect 
of which on the /3 orbit (or remnants thereof, shown as a dotted line) will depend on its strength and on the specifics of the 
superlattice (see Appendix J). The effects of bilayer splitting are not shown here for clarity, but are discussed in Appendix F. 




as one possibility, with an 
32j|), and a single variable A s 



For illustrative purposes, we adopt a model involving a single transformation of £k by Q s ,± = (tt[1 ± 25] , 7r) of 
the helical spin-density or d-density wave type such as that used in Refs 
incommensurability 6 w i « 0.06 to match incipient spin order and / or excitations 
adjusted to have the value 0.625t, as explained in the next section. The Fermi surface thus obtained comprises equal- 
sized electron (red) and hole (blue) pockets at the antinodal and nodal positions respectively, accompanied by a larger 
(green) hole pocket at the node (shown in Fig.[5]D). Here again, the line thickness is the size of the interlayer hopping. 
A significant difference in corrugation (resulting from the dispersion in k z ) is expected for pockets located at the nodal 
and antinodal locations: the obtained interlayer corrugation for each of the smaller and larger nodal pockets and the 
antinodal pocket is of order 14, 45, and 180 T respectively, being representative of general translational symm etry 
breaking order parameters yielding pockets at the antinodal and nodal regions of the Brillouin zone 12, 26- 2^, 33. 34 1. 
The difference in corrugation thus provides an opportunity for their respective locations to be identified in quantum 
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oscillation measurements. We find that the ratio of observed corrugations for the a and 7 pocket is comparable to that 
expected for the nodal and antinodal pockets. The association of the nodes and antinodes with the expected positions 
of hole and electron pockets within a broken symmetry Fermi surface 26- 2Sj]. further identify a as a two-dimensional 
hole pocket and 7 as a corrugated electron pocket within this picture. Additionally, bilayer splitting in the case where 
bilayers are ferromagnetically coupled is expected to result in a large splitting of frequencies, as opposed to the case of 
antiferromagnetic coupling between bilayers, in which case the frequency splitting of the bilayers would remain small. 
The lack of an observed splitting in any of the a, /3, and 7 frequencies is therefore indicative of antiferromagnetic 
coupling between bilayers (see Appendix F for a discussion of bilayer splitting). The observation in inelastic neutron 
scattering experiments of low energy spin excitations and a spin resonance near (71", it) that antiferromagnetically 
couple bilayers [29| suggests the association of spin order with pocket formation. 

For a larger j3 orbit corresponding to a distinct section of Fermi surface (as opposed to a magnetic breakdown 
orbit) its character can be identified in a scenario where a and 7 correspond to equal and opposite hole and electron 
pockets at nodal and antinodal locations. Luttinger's theorem (3(| (requiring a hole concentration ofp=10±l%[3|) 
and the absence of a discernible splitting of Fp from corrugation, suggests the association of the f3 orbit with a larger 
'hole-like' pocket located near (tt/2,tt/2) (see Fig. [T^.). The entire hole doping concentration would then be provided 
by the /? pocket, which acts as a reservoir of charge carriers, while the a and 7 pockets almost cancel out in carrier 
concentration. 



Electronic origin of the observed pockets 

In the simplest case, our observations indicate a broken translational symmetry order parameter that creates 
differently corrugated Fermi surface pockets at different locations in the Brillouin zone - Fermi surface models that 
yield pockets only at the nodal locations would face a considerable challenge in explaining our observations. 

Spin ordering involving a broken translational symmetry has been reported in YBa2Cu306+x samples with oxygen 
concentrations x w 0.45 [iJ H3 and « f .0 [38| in an applied magnetic field. Antiferromagnetic order has not yet 



been observed in samples with oxygen concentrations (0.49 < x < 0.8) in which quantum oscillations are observed 
12, HI HI 16 1 except when 2% of the Cu ions are replaced by spinless Zn ions in samples with x ~ 0.6 f39j. The 



strong suppression of Zeeman splitting of the Landau levels has also been reported in ref. [lj] with additional evidence 



(Appendix C) provided from the doping dependent 15j and harmonic oscillations resolved in this work, suggesting a 
translational symmetry breaking order parameter involving spin degrees of freedom or alternate means of suppression 
of the spin degrees of freedom. In addition, the implicit near-degeneracy of the bilayer split pocket frequencies (see 
Fig. [4^ and Appendices F and G) signal antiferromagnetically coupled bilayers. 

To better understand the underlying origin of the electronic structure, we map the in-plane topology of the corru- 
gated pocket. The in-plane topology of the warped 7 pocket is reflected in a variation of the in-plane calliper radius 
fen (mapped, for example in Ref. [40|, see Appendix K). In our experiment, the calliper radius is mapped by two-axes 
angle-dependent quantum oscillation measurements where 9 (the angle between H and the c-axis) is swept while 
H = H remains fixed, the sweeps being repeated for many different orientations <j) of the in-plane component of the 
magnetic field H|| = (H sm9cos(j), H sin#sin0, 0). Because of its low degree of corrugation, the method described 
here is not applicable for determination of the in-plane Fermi surface topology of the a hole pocket. The significant 
corrugation of the 7 pocket, by contrast, enables variations in fen to be mapped in Fig. [SI The quantum oscillation data 
as a function of 9 at each value of <f> is fit to a waveform composed of superimposed a and 7 oscillations as described 
in Appendix E, and the calliper radius thereby extracted at various values of <f>. Our results show a rounded-square 
in-plane topology of the 7 pocket. 

The rounded-square (i.e. nearly circular) topology of the 7-pocket (shown in Fig. Eb) is broadly consistent with an 



electron pocket predicted by numerous Fermi surface reconstruction models [12| , 26- 28, 33], although unconventional 
models - for example those involving novel quasiparticles - cannot be ruled out. The coexistence of a near circular 
electron pocket with hole pockets of comparable and/or larger size does not appear consistent with presently proposed 
'stripe' models involving multiple translations of the electronic bands £k by ordering vectors of the form Q = (ir[l ± 
2(5], 7r) where < 5 < | [33l435l ]. We consider an example of translational symmetry breaking model that could yield 
pockets of the observed size and shape to facilitate an estimation of the size of the order parameter A s that breaks 
translational symmetry. The model we consider involves a single transformation of £k by Q s .± = (7r[l ± 26], n) of the 
helical type described in the earlier section. Adjusting the value of a single parameter A s = 0.625i in such a model 
yields three pockets of size (to within 2%), carrier type, corrugation depth and shape consistent with experiment- 
see Fig. [5] and Appendix I, in conjunction with the experimentally observed constraint of hole doping p « 11.7 % 
(see Fig. [5]) and S « yg). The large ratio of A s /t indicates a substantial coupling, consistent with the near-circular 
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FIG. 6: a, Examples of angle-swept measurements made from positive to negative values of 9 on a sample of YF^CuaOe+i 
with x = 0.56 at fixed field hqH — 45 T, for different in-plane orientations (</>) of the in-plane component of the magnetic field 
H|| = (H sin 9 cos <f>, H sin 9 sin (f>, 0). Curves shown correspond to —54 < 4> < 156° (bottom to top) in steps of 15°. The gap 
in angular data around 9 — is due to a correction made for a small misalignment of the crystal, explained in Appendix A. 
b, The (^-dependence of the caliper radius k\\ of the 7 pocket determined from fits (see Appendix E) of superimposed a and 7 
oscillation waveforms (e.g. Eqn. [TJ to the data in (a)— with AF 7 ,o and F a held constant at values obtained from fits to the 
field swept data. Fit frequencies are given in Appendix E. The full 360° (f> rotation is inferred from performing a rotation over 
180° in (j) combined with fits to ±9, with only k\\ allowed to vary. The magenta line is the ^-dependence of ku expected for the 
single Q s ,± model shown in Fig. [5b. c An example fit (red line) to (f> = 51° (equivalent to 231°) data (black line), with further 
fits shown in Appendix E. 



topology of the electron pocket (the small magnitude of deviation from a circle is seen from Fig. While band 



structure calculations indicate t ~ 400 meV [18|, t : 
effective masses (see Fig. [3J lower inset), suggesting A 



100 meV provides closer correspondence with the observed 
■1 60 meV. Such a large value of A s is challenging to reconcile 
with the absence of signatures of long range order in neutron scattering experiments in YE^CuaOg+z samples of 
composition 0.49 < x ^ 0.8 (where quantum oscillations are observed). An unconventional form of broken translational 
symmetry- breaking spin order are suggested, including staggered moments present chiefly in the vortex cores [38j |. 
magnetic field-stabilised spin density wave order (37j . or spin ordering with d-wave pairing symmetry 26, 27 1; although 
unconventional models not involving translational symmetry breaking cannot be ruled out. 



III. DISCUSSION 



An open question pertains to the origin of the unexpected similarity in size of the a and 7 pockets observed in 
undcrdoped YBa2Cu306+£i;. More information on the microscopic nature of the primary instability will assist in 
addressing this question. Irrespective of their origin, the existence of compensated pockets suggests a Fermi surface 
prone to a secondary instability 0, Q at which the compensated pockets are destroyed. 
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Excitonic instability 

The possibility of an 'exciton insulator' instability (at which the compensated 7 (electron) and a (hole) pockets are 
destroyed) is raised by their similar volumes (within the experimental uncertainty, both occupying 1.91 ± 0.01 % of 
the Brillouin zone) with effective masses sufficiently large (m* ~ m* rs 1.6 m e for x — 0.56, found from the fit in the 
lower inset to Fig. EK) to cause strong attraction. Strong coupling between pairs of electrons and holes is suggested 
by a binding energy (Ryd* ~ 430 meV, derived in the Appendix B) that is more than ten times larger than their 
individual kinetic energies (-Eke ~ 39 meV). This also means that the putative exciton separation to Bohr radius ratio 
is larger than unity (i.e. r s ss 3.4 for x = 0.54, see Fig. Wp), rendering the system unstable to an excitonic insulating 
instability [41j . Due to the strongly interacting dilute Bose- Einstein condensate (BEC) limit in which this system is 
consequently positioned (i.e. r s 3> 1), the excitonic attraction is not contingent on exactly matching pocket shapes 
or perfect 'nesting.' 

Suggestive of an instability actually taking place is the occurrence of a metal-insulator quantum critical point 
(QCP) at x c — 6.46 identified in recent experiments ill 42 1. The threefold increase in effective mass on reducing 



x = 0.54 to x — 0.49 immediately preceding the QCP signals a tenfold increase in Ryd/E^ c to ~ 10 2 (see Appendix 
B) and threefold increase in r s to ~ 10 on approaching the critical doping (shown in Fig. [TJd), potentially indicating 
a further reduction in the effectiveness of screening of the Coulomb attraction, with the final trigger for an excitonic 
instability occurring at x c . Strongly correlated systems that support conditions propitious for an excitonic insulator 
are unusual 0, 0, 41 1 , given that a single carrier type rather than equal and opposite types typically results close to a 



Mott insulating state [43j. This situation appears to be reversed in YBa2Cu3 06+a;, where Fermi surface reconstruction 
at the first instability transforms the single band of carriers with a strong on-site repulsion into two bands containing 
equal and opposite pockets that are strongly attracted to each other (illustrated in Fig. [5K,b). 



Destruction of compensated pockets 



A secondary superlattice instability would be anticipated to accompany exciton condensation for x < x c [H, 
destroying the compensated a and 7 pockets (see Figs. [5J1 and e). Drawing an analogy with other underdoped 
cuprates, accompanying charge order could be expected at wavevectors Q c = 2Q S Q and / or Qj^ = (0,7r/2) [45| 
(see Appendix J). The possibility of the charge plus spin character of the second instability is suggested by the 
wavevector match Q c = 2Q S between the 'charge' wavevector Q c characterizing phonon broadenin g |4 6l| and possible 
superlattice formation [43], and the 'spin' wavevector Q s associated with short range spin order [3G| and magnetic 



excitations [311 ] . A finite charge superlattice potential A c would reconstruct the remaining /3 pocket into smaller 
pockets or open sheets [11]— the precise details depending on the choice of wavevector (see Appendix J). 

Further, the emergence of signatures of quasi-static magnetism in neutron scattering (3p|.l37j and and local moment 
behavior in muon spin rotation experiments 4t| 5(| is consistent with a strengthened A s for the region x < x c . A 
strengthened A s alone could gap the a and 7 pockets in YBaoCu^Ofi-i-.T, (see Fig.[5fc), leaving behind only the /3 hole 
pockets at the nodes - reminiscent of strong coupling models 48 1. 



Reconciliation with photoemission experiments 

Signatures from our measurements that the a pocket from which predominant oscillations arise is of the 'hole' 
carrier type located at the nodes, brings the results of quantum oscillation measurements close to reconciliation with 
the results of angle-resolved photoemission spectroscopy (ARPES) experiments. A dichotomy arises at the antinodes, 
however; while quantum oscillations indicate antinodal Fermi surface pockets measured by quantum oscillations, a 
large antinodal 'pseudogap' has been measured by ARPES measurements. An intriguing question concerns whether 
the finite spectral weight measured by ARPES to be concentrated at 'arcs' of finite extent centred at the nodal points 
(±7r/2, ±7t/2) could conceivably reflect the sole remaining parts of the j3 Fermi surface on annihilation of the equal 
and opposite a and 7 hole and electron pockets at the secondary Fermi surface instability as shown in Fig. [5] [51] — a 
Fermi surface consisting of 'arcs' is reported in materials with robust stripe order [52]. While such an interpretation is 
suggestive, an additional explanation is still required for the closing of the ARPES pseudogap in the low temperature 
and high magnetic field regime relevant to quantum oscillation measurements. 



One possibility suggested by recent models [53J, l54| for the re-emergence of spectral weight at the antinodes in 
quantum oscillation measurements is that incoherent superconducting pairs destroy the spectral weight under the 
experimental conditions (T > T c and B — 0) required for photoemission experiments, while this spectral weight is 
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recovered in the 'quantum vortex liquid' regime accessed by quantum oscillation measurements. Another contribution 
could arise from fluctuations of the excitonic insulating instability that cause loss of antinodal spectral weight over a 
range of oxygen concentrations above x c extending to higher temperatures— similar to the case of NbSe2 and TiSe 2 
above their charge ordering temperatures 57, 58|. The possible role of density wave reconstruction has also been 
suggested by theory (HBJ and by results of ARPES experiments [H^, [6(| ■ Alternatively, more exotic models involving 
novel quasiparticles visible to quantum oscillations and not ARPES have also been proposed [56 1. 



Relevance to pairing 

Finally, we consider the effect that excitations of the potential excitonic insulator instability may have on the 
quasiparticle pairing mechanism. The unusual occurrence of a Fermi surface ideally predisposed to an excitonic 
insulator instability in the same material where high T c superconductivity occurs could either suggest an link between 
the phenomena, or a coincidence. Similar compensated corrugated cylindrical Fermi surface pockets are found in 
other families of unconventional superconductor that also exhibit unexpectedly high transition temperatures— such 
as the Fe pnictides [61] and Pu-based superconductors [62]. 

Theoretical proposals made prior to the advent of high temperature superconductivity considered conditions under 
which an excitonic mechanism may potentially enhance T c Q . Whether or not such a cooperative mechanism is in play 
in the d-wave cuprate superconductors, however, is contingent on a better understanding of the relative importance 
of spin and charge fluctuations in contributing to high T c in these materials. 



APPENDIX 



A. Experimental details 



Quantum oscillations are measured using the contactless conductivity method 63[ in two pulsed magnets reaching 
65 T and 85 T, and in a DC magnet reaching 45 T. Plate-like single crystals of YBa2Cu3C>6+x with in-plane 
dimensions 0.3 x 0.8 mm 2 are attached to the face of a flat 5 turn coil (compensated in the case of pulsed field 
experiments) that forms part of a proximity detector circuit resonating at ~ 22 MHz [gij]. A change in the sample 
skin depth (or complex penetration depth in the superconducting state 65]) leads to a change in the inductance of 
the coil, which in turn alters the resonance frequency of the circuit. Quantum oscillations are observed in the in-plane 
resistivity and hence the skin depth). 

To minimize the effects of flux dissipation heating in pulsed magnetic fields, the sample is immersed in liquid 
4 He throughout the experiment 0, [67} , and only the oscillations observed during the falling field are considered for 
detailed analysis. On taking such precautions, effective mass estimates are obtained that are comparable to those 
obtained in samples of the same oxygen concentration in static magnetic fields [12] • 

The fixed H (= 45 T) angle swept 9 measurements performed at different in-plane angles <fi shown in Fig. [6] are 
obtained using a dual-axis rotator. Due to a small misalignment of the sample by 9q « 7° at cj>o w 51°, a minor 

correction cos factual = cos (^uncorrected x \ cos 2 9q + sin 2 9q sin 2 ((/) — 0q) is made to obtain the correct 9. 



B. Estimates of the exciton binding energy 

We estimate the effective 'Rydberg' for exciton binding to be 0, l4l| 

Ryd* = -^^Ryd (7) 

k, 430 meV for YP^CuaOg+z samples of oxygen concentration x = 0.54 and 0.56 [7J, where \i sa m*/2 sa 0.8 m e 
is the reduced bi-exciton mass, Ryd ss 13.6 eV is the hydrogen Rydberg and e ~ 5 is the relative permittivity for 
YBa2Cu306+a; [68]. The kinetic energy is estimated using 



(8) 
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« 39 meV for F — 535T and m* = 1.6 m e , yielding a ratio Ryd* /E^ « 11. The effective Bohr radius of the exciton 
is estimated to be 

a = e — ao (9) 
A 1 

rs 3.3 A, where ao « 0.53 A is the hydrogen Bohr radius. We estimate the exciton density as n = 
Ak/2TT 2 w 2.6 x 10 17 m 2 , where Afc = 2neF/h rs 5.1 x 10 17 m _1 , yielding an effective exciton separation of 



(10) 



1 / 1 
a* V ftn 
from which r s w 3.4. 

In estimating i?yd* and r s , we have neglected the effects of screening, which could become significant should a large 
(3 pocket survive the formation of a superstructure accompanying exciton condensation. On the other hand, super- 
structure formation (see following section) would likely cause significant reconstruction of the f3 orbit— the outcome 
being strongly dependent on the associated ordering vectors and nature of the broken symmetry (see Appendix J). 
Screening will also become ineffective if to* globally diverges at the quantum critical point at x c 0.46, at which 
point both Ryd* and r s would diverge while E^ c collapses to zero. 

When the quasiparticle effective mass changes (as reported in YBa2Cu306+a; as a function of doping 15]), then we 
can expect the following dependences (shown in Fig. Q] of the main text): 

Ryd* 



Eke 



oc (to*) and r s oc to* (11) 



indicating that tuning of the electron correlations is expected to have dramatic consequences for the likelihood of 
exciton pairing and the applicability of the strong coupling BEC limit [l], |4l[ . 



C. Extended analysis of spin splitting 

For a finite value of 7 (which is the effective moment of the spin doublet (i.e.) the Wilson ratio in the param- 
agnetic phase), the spin damping factor R s — cosOjrym* /m e ) has an oscillatory form. In the case of underdoped 



YBa2Cu306+x, however, it has been shown in ref. [14| and in greater detail here, that the value of 7 <C 1 for the a 
pocket corresponding to the spectrally dominant oscillations, leading to a value R s w 1. The effect of an oscillatory R s 
would be for the quantum oscillation amplitude to vanish at 'spin zero' values of ym*/m e = (2n + 1) /2, on either side 
of which the phase of the quantum oscillations changes by 7r, resulting in their inversion. We test for an oscillatory 



R s by examining the phase of measured quantum oscillations for a range of to* tuned by angular rotation [14[ and by 
changing the doping [15j. Figure [7^,, b shows the phase of the fundamental oscillations to be invariant as a function of 
to* from 1.7 to 3.0 to c on angular rotation, and from to* from 1.6 to 4.5 m e on varying the doping. Figure[7fc shows 
the phase of the harmonic (2F) oscillations to be invariant as a function of to* from 3.4 to 4.4 to c on angular rotation. 
The measured quantum oscillations dominated by the a frequency are seen to be phase invariant for a wide range in 
to*. From figure [TJi, a finite 7 would lead to multiple spin zeros over the measured range in to*, the absence of which 
signals a value of 7 <C 1 for the dominant a frequency. Consequently, R s rs 1 does not appear in Eqn. [TJ 



D. Extended analysis of swept magnetic field data 

During Fourier analysis, weaker spectral features can be affected by interference effects between oscillations from 
a spectrally dominant series of oscillations (nearby in frequency) and the window function (a Hann window in the 
current study) used to modulate the data. Such interference between the a oscillation and the window function affects 
the resolved 7 frequencies in the present study. One way to mitigate interference is to subtract fits to a dominant well 
characterized sequence of oscillations (i.e. the a oscillations) prior to Fourier analysis. This procedure is justified in 
the present case owing to the resolved separation between the a and 7 frequencies in our Fourier analysis (see Fig. 31 
Fig. [3l and Fig. with the degree of separation growing proportionately larger for the harmonics as expected. 

In Fig. [8] we show Fourier transforms of the measured oscillations in samples with x — 0.54 before (i.e. Fig. |4]d of 
the main text) and after (i.e. Fig. 0fc of the main text) subtracting the fit of the a oscillations (Eqn. Q] in the main 
text), a subset of which are shown in Figs.[3jD and c of the main text. The 7 frequencies -F-y.ncck and F 7i b e iiy in Fig.Hk 
of the main text are determined with greater accuracy by performing Fourier transformation after subtracting the a 
oscillations. 
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FIG. 7: a Fundamental quantum oscillations measured for a range of effective masses tuned by doping, taken from ref. [151 ]. 
Here the quantum oscillation amplitude has been rescaled by an exponential factor for purposes of comparison, b Fundamental 
quantum oscillations measured for different effective masses tuned by angular rotation, taken from ref. [TJ]. Each set of 
quantum oscillations is labelled by the tuning parameter (angle, doping) and the corresponding effective mass marked by 
arrows. The gray vertical lines label the expected peak positions for the fundamental a quantum oscillations in b and c. c 
Harmonic quantum oscillations measured for a range of effective masses tuned by angular rotation, obtained by subtracting 
the fundamental alpha and gamma oscillations from the overall signal measured in the current work (Fig. |4p). Vertical cyan 
lines label the expected peak positions for the harmonic (2F) a quantum oscillations, a, b, and c reveal an absence of phase 
inversion of the measured quantum oscillations over a broad range in effective mass, d The effective mass locations of the spin 
zeros expected for different values of the spin doublet moment 7. The value of 7 (i.e. the Wilson ratio) in the paramagnetic 
phase is marked by the horizontal blue line. The yellow shading represents the range of effective mass tuned by the doping 
in a, the cyan shading the range of effective mass tuned by the angle in b, and the pink shading the range of effective mass 
covered by the angle-tuned harmonics in c. Given the absence of phase inversion (which would locate a spin zero) in quantum 
oscillations over the measured range of effective mass, 7 <C 1 is indicated for the spectrally dominant a pocket. 



E. Extended analysis of swept angle data 

When a Fermi surface has an interlayer corrugation, as we identify for the 7 pocket (depicted in Fig. [1^), angle- 
dependent measurements provide a means for mapping its in-plane topology [20j. If the degree of corrugation is 
sufficiently deep to yield separately resolved neck and belly frequencies in Fourier transforms [23[ , quantum oscillations 
then provide an accurate means for mapping the in-plane topology. The in-plane topology can be efficiently mapped 
out by sweeping 9 in a constant magnetic field for many different in-plane orientations <f> of the magnetic field 
H|| = (H smO cos(/>, H sin0 sin 0, 0)— requiring the use of a dual-axis rotator (for other examples, see Appendix K). 
In the present experiments on YBa2Cu306+ K with x — 0.56 and 0.54, the 7 electron pocket is found to be deeply 
corrugated yielding spectral features between 10 and 40 % in amplitude of the spectral weight associated with the a 
pocket. A minimal model for fitting is given by 

A/ = Af at0 cos — exp + A/ 7 , cos — - exp 

\HoH cost) / \ fi H cos u J \[ioH cost) J \ hqHcosOJ 

2^AF 7|0 J (fc||Ctanfl) \ 

HoHcosO J [ ' 

(similar to Eqn. [IJ but without including a corrugation of a) that takes into consideration the superposition of a and 
7 oscillation waveforms in Figs.lHk. an d E^). Here A/ Qi o and A/ 7i o are amplitude prefactors, T 7 is a damping factor 
for the 7 oscillations, while the remaining parameters are defined in the main text. 
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FIG. 8: a, b, c and d, Fourier transforms of the oscillations measured in YBa2Cua06+a; with x = 0.54 at different angles 9 
before (mauve line) and after (magenta line with grey shading) subtracting the a component of a fit to Eqn. Q] of the main 
text— a subset of which are shown in Fig. [3]d and c of the main text. The resolved 7 frequencies are used to construct Fig. [4^, 
of the main text. 




FIG. 9: a Swept 9 oscillations measured in a sample of YF^CusOe+i with x = 0.56 at constant field H = 45 T and different 
fixed values of at T ~ 1.5 K fit to Eqn. [12] (black lines) as described in Fig. [6] The gaps in the data result from a correction 
for a sample misalignment of ta 7° (see experimental details), b Fourier transforms in 1/fioH cos 9 of the oscillations shown in 
a, which are dominated by the a pocket. As expected for a nearly ideal two-dimensional Fermi surface section, the dominant 
a frequency obtained is the same as that from fixed 9 swept magnetic field experiments (i.e. Fig. [6^0 at all angles of rotation. 
The curves correspond from bottom to top to —54 < <j> < 156° in 15° steps. 



Changes in the caliper radius of the 7 pocket fcii are assumed to be the dominant factor responsible for the subtle 
</>- dependence of the waveform in Fig. |6^i of the main text and Fig. [9^i. Fourier analysis in Fig. [9}d of the (9-swept 
oscillations in I/hqHcosO yields a leading frequency F a w 539 ± 5 T similar to that obtained on sweeping the field 
at fixed angle, confirming the l/cos# angular dependence of the dominant a oscillations. We note that an additional 
angular damping mainly affecting the 7 pocket is observed on comparison of the oscillations measured as a function 
of angle with the oscillations measured as a function of magnetic field — suggesting an interlayer decoherence factor 
that additionally suppresses the oscillation amplitude with increased angle. 

Owing to the subtle variation of the waveform with cj>, parameters Af a ^, A/ 7i o, F at o, Kyfl and T 7 are obtained 
from fits of Eqn. [12] to the angular data assuming a circular cross-section (i.e k\\ = fcp, where fcp = y / '2eF~ ( ^/h is the 
mean Fermi radius). The fits yield frequency parameters F Qi o = 536 ± 5 T and F 7) o = 535 ± 5 T, closely matching 
those determined in magnetic field sweeps. These parameters are then held constant for all 4>, and a single parameter 
fit is performed with only k\\ being allowed to vary to refine the fits for each <j>. The resulting ^-dependence of feii — 
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shown in Fig. [6p of the main text— yields a rounded-square Fermi surface cross-section for the 7 electron pocket, 
similar to that calculated using the helical spin-density wave model. We assume that the in-plane pocket topology 
accessed is an average between the bilayer split pocket topologies (such as shown in Fig.[TU]). In Fig. [5^,, we show the 
4> dependence of the actual fits to Eqn. [FZ[ 



F. Effect of bilayer splitting 



Bilayer splitting in YE^CuaOg+x arises from the bilayer crystalline structure, yielding bonding and antibonding 
bands that lead to two large hole sections of different sizes in the paramagnetic phase (see Fig. [TOk ) fljj . Here the 
width of the line represents the difference in area between bilayer-split bands, with '+' and ' — ' representing the sign 
of this difference. 

If ordering occurs in which there is antiferromagnetic coupling between bonding and antibonding bands, the Hamil- 
tonian factorizes into 



-a 

-k+Q 



and 



A F h 



(13) 



The resulting Fermi surface shown in Fig. 110b has two different variants— the difference between them being more 
pronounced for the electron pocket, which is stretched along each of the k x and k y directions. The difference in area 
between each is represented by the thickness of the line, with the sign of the difference being positive or negative as 
indicated by '+' or ' — .' Since the lines bounding each of the pockets contain both positive and negative segments 
that almost cancel, the pocket areas of the each bilayer-split variants of the reconstructed Fermi surface are only 
weakly split— i.e. the degeneracy of the pockets is not significantly lifted by bilayer splitting when the bilayers 
are antiferromagnetically coupled. Consequently, different frequencies would not be expected to arise in quantum 
oscillation experiments [26l . l27j | from antiferromagnetically coupled bilayer split pockets. 



(a) (7i,0) 




(b) (7i,0) 

l T 




FIG. 10: Calculated Fermi surfaces in which the effects of bilayer splitting are included as a perturbation in r (here the effects 
of corrugation are not shown for clarity, but are the same as in Fig. [5]). The difference in area between different bilayer-split 
variants of the Fermi surface is represented by the thickness of the line, where the colours (and '+' and ' — ' labels) indicate 
the sign of this difference, a shows the translation of the unreconstructed Fermi surface where the coupling between bonding 
and antibonding bands is represented by different signs for the difference in area (black for ' — ' and gray for ' + ') between £k 
and £k+Q- b shows the reconstructed Fermi surface where the difference in area (red for ' + ' and pink for ' — ' for the electron 
pocket, and blue for ' + ' and cyan for ' — ' for the hole pocket), can be seen to cancel for each pocket. 



Were the bilayers instead ferromagnetically coupled, then like bands would be coupled giving rise to significant 
differences in pocket areas between each of the bilayer-split Fermi surfaces. The size of this difference is expected to 
significantly exceed the depth of corrugation (l8j and so would be observable as distinct frequencies in the quantum 
oscillation data, observable as well separated peaks in the Fourier transform. 



G. Inconsistency of scenario where observed multiple frequencies arise from bilayer-split pockets 

Here we show that a scenario where the observed multiple frequencies arise from two bilayer split pockets such as 
invoked in Ref. (l3| is inconsistent with the experimental data. Fig. [11] shows the least square error in an unconstrained 
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fit of the measured oscillations to Eqn. [IJ The Bessel term Jo(2TrAFi/B) that modulates the quantum oscillation 
amplitude in Eqn.Q] 21] arises from phase smearing due to the Onsager phase of a quasi- two-dimensional Fermi surface 
with depth of corrugation AFi that varies sinusoidally with k z . In the limit 2nA.Fi/B >> 1/4, the waveform from each 
section i is a superposition of extremal neck and belly frequencies Fi neck,beily — Fi -F Ai^. For corrugated cylindrical 
pockets with a small degree of warping, such as the a pocket, the difference in Onsager phase (2tt(2AF)/ ^i n H) 
between the neck and belly extremal orbits at (j,qH = 65T is small enough (« 2 radians for the a pocket) to render 
inadequate the parabolic approximation assumed in the Lifshitz Kosevich expression (22| . Instead, the corrugated 
cylinder topology is best represented by Eqn. [T] in which the amplitude of the neck and belly oscillations is constrained 
to be identical (i.e. Cj), and their relative phase difference constrained to be ir/2 2l|. 

All parameters have been allowed to vary for different fixed ratios of AF 7i o and AF a ^, and the resulting least 
squares error of the best fit obtained as a function of the corrugation ratio between the two pockets. It is evident 
that the minimum error appears for a significantly different corrugation of the a and 7 sections, with a rapidly 
growing error with increasingly similar corrugations of the a and 7 Fermi surface sections. Very similar parameters 
are also obtained on fitting the quantum oscillations in Ref. [HI to Eqn. Q] (where a sine term replaces the cosine 
term for magnetic oscillations), indicating consistency with two very differently warped pockets. Considering the 
scenario where bilayer split pockets are invoked to explain the multiple low frequencies, the warping of each of the 
two a and 7 sections would be required to be the same (we assume that the non-conducting chains in underdoped 
YBa 2 Cu 3 06+x do not hybridise with bonding and anti-bonding states to create differently warped bilayer split sections 
as in optimally doped YBa2Cu307 flij]). The least squares error is near maximum in this case, as demonstrated by 
the best fit for AF 7i o/ AF at o « 1 shown in Fig.QTJ This disagreement in fit indicates that the multiple observed low 
frequencies are not explained by similarly warped pockets, as would be the case for bilayer-split pockets originating 
from a ferromagnetic bilayer coupling. 

The bilayer splitting which must occur for YBa2Cu306+x is not experimentally resolvable as an additional set of 
frequencies in the Fourier transform. Indications therefore are that the areas of the bilayer split pockets remain almost 
degenerate, as would occur for the antiferromagnetic coupling of bonding and antibonding states, consistent with the 
low energy spin excitations observed in inelastic neutron scattering experiments 29] . 
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FIG. 11: Error analysis of unconstrained oscillation fit to two Fermi surface sections, (a) Least squares error of the unconstrained 
fit of Eqn. \T\ to the measured oscillations of the sample with oxygen concentration x = 0.56 as a function of different ratios 
between the corrugation of the a and 7 pockets is shown on the left hand axis. The corresponding difference between the fit 
corrugations of the a and 7 pockets is shown on the right hand axis. A clear minimum is observed for AKy t o/AF a ,o ~ 6.2, 
yielding the best fit parameters - while the error is seen to rapidly increase if the corrugations of the two pockets are similar, 
(b) Example best fit (red line) to the measured oscillations for the scenario where AF lt o ~ AF a .o. The fit does not yield good 
agreement, demonstrating the need for two low frequency Fermi surface sections of significantly different warpings to explain 
the measured oscillations. 



It can readily be seen why the scenario of two equally warped Fermi surface sections does not explain the multiple 
low frequency oscillations. We define two pockets (labeled T' and '2') with neck and belly frequencies -Fi, n eck and 
Fi, belly and F2,neck and F^bciiy in terms of two median frequencies F\ and F2 with a depth of corrugation AF: i.e. 

^1, neck/belly = Fl ± AF , i^.ncck/bclly = i*2 ± AF. (14) 

where AF would have the same magnitude and sign for the two pockets in a ferromagnetically-coupled bilayer splitting 
model. 
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Quantum oscillations result from a superposition of two oscillatory terms 



A/ oc 



2ttFi 
[i H cos 9 



Jo 



2ttAF 
Mo-ff cos 9 



2ttF 2 



•Jo 



/ 2k&F 



/i a Hcos9j \HqH cos 



for the two pockets, yielding the four frequencies given above. On combining them, we obtain 



A/ oc cos 



fio H cos 9 



tt(Fi - F 2 )\ / 2ttA^ 



Jo 



/iqHcos9 J \fJ.oH cos 9 



(15) 



(16) 



Two terms in Eqn. \W\ will give rise to nodes: the Bessel term containing AF will give rise to nodes at the set of 
magnetic fields given by Eqn. |3] (eqn. [3] in the main text), while the cosine term containing (Fi — F%)/2 will give rise 
to nodes 



fi Q H m cost 



\Fi — Fi 
2m + 1 



(17) 



where m is an integer. Since the node terms are multiplicative, both sets of nodes would be observed giving rise to a 
suppression of the total waveform (as shown in the attempted fit in Fig. ITTb ) at the field values given by both Eqns. [3] 
and [T7J The absence of such nodes in the full experimental waveform (Fig. 2a, 0t> of the main text) indicates that 
the multiple observed low frequencies are not a consequence of bilayer-split pockets. 



H. Additional swept angle data 

In Fig. [12] we show oscillation in data measured by sweeping 9 for constant H at many different values of </> obtained 
for a single crystal sample of YBa2Cu3 06+x with x — 0.54. The weak </>- dependence reflects the behavior observed 
for x — 0.56. 
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FIG. 12: Swept 9 oscillations measured in a sample of YBa2Cu306+a; with x = 0.54 at constant field H = 45 T and different 
fixed values of d> at T ~ 1.5 K. 



I. Bandstructure details 



Antiferromagnetic coupling between bilayers couples bonding to antibonding bands, yielding two slightly different 
variants of the Fermi surfaces shown in Figs.[5]D and [5J; of the main text with very similar areas. 

For samples with oxygen concentration close to x w 0.5, ortho-II ordering has the potential to modify the Fermi 
surface, depending on the strength of the relevant potential Aj t 2ol 2 
An alone can yield pockets comparable to those measured [2_ 

ortho-II makes little difference to the observed frequencies [27J (see Fig. [T5|) . In Fig. [T3] the ortho-II potential can 
been shown to provide improved agreement between the calculated Fermi surface and experiment. 



For commensurate spin ordering (i.e. 6 = 0), 
For incommensurate ordering where 6 > jg, however, 
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FIG. 13: The left-hand panel shows the Fermi surface calculated for a helical spin-density wave model. Using a hole doping 
p — 11.7 % and 8 ~ ~ 0.06 as constraints, adjustment of a single parameter A s yields a Fermi surface with frequencies 
comparable to those observed— i.e. F a fi = 575 T, Fp t o = 1583 T, and F 7 ,o = 517 T. The largest discrepancy for the case of the 
/3 pocket corresponds to 0.2 % of the paramagnetic Brillouin zone area. The helical spin-density wave model used is described 
in Ref. fl2l ]. On including an orthoTI potential (right hand panel), using the formalism described in Ref. [2q|. pockets with 
frequencies— F a ,o = 517 T, Fp,o = 1641 T, and F 7i o = 515 T— closer to those observed experimentally are obtained in the 
right-hand panel. The largest discrepancy for any one pocket now corresponds to 0.07 % of the Brillouin zone. 
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FIG. 14: The actual Brillouin zone for Fermi surface reconstruction by Q a = (ir[l ± 2<5],7r), assuming 8 = y^, illustrating the 
possible superlattice modulation vectors Q c and Qx accompanying an excitonic insulator instability. 



J. Possible ordering resulting from an excitonic insulator instability 



While the extended zone representations in Figs.[5jD and c of the main text provide a convenient means of visualizing 
the three distinct pockets, the actual magnetic Brillouin zone depicted in Fig. IT4l (obtained on repeatedly folding that 
in Fig. [5}d of the main text) must be considered in order to identify possible ordering vectors associated with the 
excitonic insulator instability. Here, we make the simplifying assumption that S = A ~ 0.06. 

One possible candidate ordering vector is Q c = 2Q S = (±45ir, 0), corresponding to the charge modulation that can 
accompany a spin-density wave should it become collinear or 'stripe-like' (as opposed to helical). In this case Q c does 
not introduce a new periodicity in Fig. [14j but rather introduces additional couplings between £k+nQ a and £k+(n±2)Q s 
that are not present in a purely helical spin-density wave. Such additional couplings will cause gaps to open around 
the edges of the magnetic Brillouin zone in Fig. [14j leading to the destruction of closed pockets and the creation of 
open Fermi surface sheets in a manner analogous to those obtained Ref. [HI (where S = h was considered). The 



observation of phonon broadening effects [46] together with the possible observation of Bragg peaks in Ref. 47] could 
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be consistent with charge ordering of this type. 

An alternative possibility is a form of broken translational symmetry with a characteristic vector w (0, tt/2) 
that maps the electron and hole pockets onto each other as in a conventional Fermi surface nes ting scenario. In this 
case, Qj^ has a similar periodicity to the purely charge wavevectors identified in stripe systems [331 ] and/or scanning 
tunneling microscopy experiments 45|. While such forms of order are reported in various families of high T c cuprates, 
they have not been observed in YBa2Cu306+i. 
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FIG. 15: a, Experimentally measured angular dependence of the dHvA frequencies in LaRhlns. b Band structure calculations 
of the Fermi surface topology of the a-pockets in LaRhIn5 



K. Rotation studies on other corrugated Fermi surface materials 

The rare-eath paramagnet LaRhlns and the organic superconductor /3-(BEDT-TTF)2lBr2 are two examples 
of strongly correlated materials for which a similar degree of corrugation enables angle-dependent quantum oscillation 
experiments to access information on the in-plane topology. 




Angle (deg) 

FIG. 16: a, In-plane Fermi surface cross-section of /3-(BEDT-TTF)2lBr2 from Ref. [2(j. b, ^-dependence of the difference 
frequency Af\,3 = Fbelly — -fncck (also from Ref. [20j]) for two different <j> angles, which 90° apart according to (a). 

LaRhlns has a section of Fermi surface (labelled a in Fig. 1 15[) with a proportionately similar degree of corrugation 
to the 7 pocket in YBa2Cu306+z and a round-square cross-section de Haas- van Alphen experiments can resolve 
neck and belly frequencies F a i^ and F a ^ respectively, which become degenerate when H is rotated by an angle 9 
away from the [0,0,1] axis, corresponding to the first Bessel zero in Fcos9 ps F ± AFJ (k\\c' tan0). Owing to 
the rounded-square cross-section, the caliper radius is larger along < 1,1,0 > than along < 1,0,0 > causing the 
degeneracy to occur at a smaller angle 9 when the field is rotated towards [1,1,0] than when it is rotated towards 
[1,0,0]. 

According to Ref. [20j, /3-(BEDT-TTF)2lBr2 has a warped Fermi surface (see Fig. ITBk) with an elliptical in-plane 
topology and triclinic crystal structure, causing the ^-dependences of the difference frequency AF1.3 = Fbciiy — -Fb ec k 
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(see Fig. [TBb) to be different for the two different <j> angles considered (which are 90° apart). The existence of 
only a single section of Fermi surface /3-(BEDT-TTF)2iBr2 and the very low transition temperature of this organic 
superconductor enables multiple nodes to be directly observed. 
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